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The multivariate extremal index function relates the asymptotic distribution of the vector of 
pointwise maxima of a multivariate stationary sequence to that of the independent sequence 
from the same stationary distribution. It also measures the degree of clustering of extremes in 
the multivariate process. In this paper, we construct nonparametric estimators of this function 
and prove their asymptotic normality under long-range dependence and moment conditions. 
The results are illustrated by means of a simulation study. 
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1. Introduction 

The motivation for this paper comes from an empirical observation that time series from 
hydrology, meteorology, environmental sciences, finance, etc. are heavy-tailed and clus- 
tered when extremal events occur. In particular, it has been recognized in recent decades 
that the model of independent and identically distributed (i.i.d.) Gaussian random vari- 
ables is inappropriate for modeling extreme returns of risky assets that are observed 
during a financial crisis. It is important for risk managers to understand the relative be- 
havior of the various financial risks to which their institutions are exposed in the event of 
large losses because they have to anticipate the diversification opportunities so that the 
risks can be balanced by comovements (between risks) or reversal movements in short 
time intervals (within risks). 

Although there are well-developed statistical approaches to characterize the cross- 
sectional dependence structure of extreme returns of risky assets (see, e.g., [14, 20, 23, 
34] and the references therein), problems concerning the estimation of their temporal 
dependence structure have not received much attention. A notable exception is [46], which 
proposes a specific class of max-stable processes to model simultaneous dependencies 
between and within financial time series. However, this ad hoc class of processes is not 
necessarily suitable for any multivariate time series. The multivariate extremal index 
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function, introduced by Nandagopalan [27, 28], is a quantity which allows one to relate 
the asymptotic distribution of the vector of pointwise maxima of a stationary sequence to 
that of the independent sequence from the same stationary distribution. It also measures 
the degree of clustering of extremes in the multivariate process since it is equal to the 
reciprocal of the mean number of clustered extremal events. Therefore, it is a specific 
measure of the temporal dependence structure of the extreme values of the process. 

It is the aim of this paper to present a general theory for the inference of this function. 
We extend the block declustering approach introduced in [37] to the case of multivariate 
stationary processes: we construct pointwise estimators and study their asymptotic prop- 
erties. Three assumptions are made: (i) there exist moment restrictions on the amount 
of clustering of extremes; (ii) the number of two-level exceedances converges weakly - 
an assumption which will guarantee the existence of the asymptotic variancc-covariance 
matrix of the estimators; (iii) a mixing condition weaker than strong mixing is supposed 
to hold. Under these assumptions, we prove the asymptotic normality of our estimators. 

More formally, let (Xi = (-Xj i, . . . , Xi y d))i>i be a strictly stationary sequence with sta- 
tionary distribution function F{x) = P{Xi : i < x%, i = 1, . . . , d), x = (xi, . . . , Xd) G M. d , and 
univariate marginal distributions Fi(x) = P(Xij < x), i = 1, . . . , d. We assume that there 
exists a family of normalizing sequences in M d , (m„(t) = (u n! i(Ti), . . . ,u ny d(Td))) n >i, 
t = (n, . . . , Td) G (0, oo) d , such that 

lim n(l — Fi{u n _i{T))) = r for r > 0, i = 1, . . . , d, (1.1) 

n — >oo 

and, for some function H : (0, oo) d i— > [0, 1] , 

lim n(l - F(u n (r))) = - In H(t), for r G (0,oo) d . (1.2) 

n — >oc 

A necessary and sufficient condition for the existence of a sequence (u Ili? :(r)) rl >i which 
satisfies (1.1) is that lim x -, Xf i F(x)/Fi(x— ) = 1, where Xfj = sup{w : Fi(u) < 1} and 
F = 1 — Fi (see Theorem 1.7.13 in [22]). A natural choice for u nt i(r) is then given by 
F?~ (1 — i~/n), t G [0,n), where F*~ is the generalised inverse of Fi, that is, F*^(t) — 
inf{cc G R : Fi(x) > r}. This assumption is weaker than assuming that Fi is in the domain 
of attraction of an extreme value distribution since the normalization is linear in this 
case. However, the function G defined by G(r) = ^(r-f 1 , . . . jT^ 1 ) for -r G (0,oo) d must 
be a multivariate extreme value distribution regardless of whether the normalization is 
linear (see [33], Proposition 2.1). In particular, G is a continuous distribution function 
with unit Frechet margins. It is noteworthy that — In If is a homogeneous function of 
degree 1, that is, — lniJ(cr) = — clnij(r) for all c > and r G (0, oo) d . This function is 
sometimes called the stable tail dependence function of F . 

Let M n> i = max(A"i i i, . . . ,X n ^) be the maximum of the ith component and introduce 
the vector of pointwise maxima M n = (Af„.i, . . . , M„^)- If (X„) n >i is a sequence of 
independent and identically distributed (i.i.d.) vectors of random variables (r.v.s), then 
(1.2) is equivalent to 

P{M n < u n (r)) = P(M n>i < iLn,i(n)>i= 1, • --,d) -> H(t), as n—> oo. 
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This convergence can be extended to stationary sequences by assuming the long-range 
dependence Z)(M n (r))-condition introduced in [19], which is a natural multivariate ver- 
sion of the well-known univariate Z?(u„(r))-condition (sec, e.g., [22], page 53). Let 
§+ = {t e (0, oo)'' : \t\ = 1}, where |r| 2 = £? =1 if. If D(u„(t)) holds for each r € (0, oo) rf 
and P(M n < m„(to)) converges as n — > oo for each To in §+, then there exists a function 
9 : (0, oo) d i > [0, 1] such that (see Proposition 2.3 in [33]) 

lim P(M n <u n M) =H{t) 6[t) for r G (0, oo) d . (1.3) 

n — >oc 

Let if (r) = F(r) fl W. The function G defined by G(r) = F(rf x , . . .,-r^ 1 ) is also a mul- 
tivariate extreme value distribution and has Frechet marginals. The function 

is referred to as the multivariate extremal index function of (-X'„)„> 1 . 

The estimation of the function — In if for sequences of i.i.d. vectors of r.v.s has been 
thoroughly investigated (see, e.g., [9, 10, 11, 16, 18]). But, extensions to stationary se- 
quences are still at an early stage (see an example in [42]). The estimation of the mul- 
tivariate extremal index function has been little investigated. Recently, some pointwise 
estimators have been proposed, but their asymptotic properties have not been studied. 
In [41] , Smith and Weissman introduce the class of multivariate maxima of moving max- 
ima (Mi) processes and establish that the multivariate extremal index function of a very 
wide class of processes may be approximated arbitrarily closely by one from a M4 pro- 
cess. Nevertheless, the estimation of 8 via a M4 process is practically infeasible since 
it necessitates the estimation of an infinite number of parameters, except if additional 
approximations are made. Smith and Weissman also give a key characterization of 6(t) 
as the univariate extremal index of a sequence depending on the standardized Frechet 
components (see Proposition 2.1 in [41] and Proposition 2.1 below). It follows that if one 
can transform the data to have unit Frechet components, then 9(t) may be estimated 
by univariate methods such as those proposed in [13, 40, 45], or [37]. To evaluate the 
accuracy of this approach, a simulation study is conducted in [12] with Ferro and Segers 
estimators (see [13]). 

In this paper, we introduce two new nonparametric estimators of the multivariate 
extremal index function. Its original contribution is to study the asymptotic properties 
of these estimators. The paper is organized as follows. In Section 2, we discuss two 
characterizations of the multivariate extremal index function and present some of its 
properties. In Section 3, we explain how we construct the estimators. Note that they 
are based on a block declustering scheme and are only determined by the block length, 
as in [37]. In Section 4, we present and discuss technical conditions. We then give the 
asymptotic distributions of the estimators. In Section 5, we investigate their finite-sample 
behaviors on simulated data. The proof of the asymptotic normality of the estimators is 
found in Section 6. Section 7 concludes. 
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2. The multivariate extremal index 

The multivariate extremal index function defined by (1.4) may also be characterized by 
the asymptotic distribution of the following point process of exceedances: 



N i T H B ) — 1 {i/neB,X ; ^u„(-r)}> 
1=1 

where B is a Borel set included in (0,1] and {X; ^ u n (r)} = U»=i{-^M > u n,i( T i)}- 
Contrary to the univariate case, there are several ways to define a multivariate threshold 
exceedance. Here, we define an exceedance as the event that one of the components of Xi 
exceeds its associated threshold. Suppose that (1.3) holds and (X n ) n >i satisfies the long- 
range dependence A(it„(r))-condition introduced in [28] (which is a little stronger than 
the Z?(M„(r))-condition). A necessary and sufficient condition for the weak convergence 
of Nt\-) is then the convergence of N^' ((Q;q n /ri\) to adiscrete distribution, 7r^ T \ given 
that there is at least one exceedance, that is, 

lim P(N^((0;q n /n]) = k\N^((0;q n /n})>0)=^\k), 

n — >oo 

where (q n ) is a A(u n (r))-separating sequence (see Section 4). ir 1 -^ is referred to as 
the cluster-size distribution. Under these assumptions, the point process converges to 
an homogeneous compound Poisson process, N^ T \ with intensity — 9{t) lniJ(r) and 
limiting compound distribution . It may be stressed that, under some mild additional 
assumptions, we have the following characterization of the multivariate extremal index 
function (see [19, 24]): 



9{t)= X> W (*0 
\fc=i 

that is, 9 is equal to the reciprocal of the limiting mean number of exceedances in a 
cluster. 

As mentioned in the Introduction, Smith and Weissman give an alternative characteri- 
zation of the multivariate extremal index function in [41] . They first propose to standard- 
ize the margins to the unit Frechet distribution and then to express 9(t) as the univariate 
extremal index of the constructed sequence as a linear combination of the standardized 
components. In this paper, we decide to standardize to the unit Pareto distribution, as 
in Section 10.5.2 of [4] (see property (v)). 

Proposition 2.1. Let r <G (0,oo) d \{0} and assume that lim rwoo P(7W„ < u n (r)) = 
H{t) 9 ( t \ Define the associated univariate stationary sequence by 
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Z, = max TiYn, l>l, 

1 i=l,...,d 

where Yij = (1 — (Xij))^ 1 and (aj) = P(Xu < x). 0(t) is then the univariate 
extremal index of the sequence (Zj )/>i, that is, it satisfies, for x> 0, 

lim nP[z\ T) >v£ ] {x)) = x and lira. P[ max Z, M < v<£ ] (x) ) = c" 0(x)5< , 

n — >oc n — >oo \ /— l,...,n 7 

where v^\x) = >r _1 (— lnff(T))n. 

The proof of Proposition 2.1 is postponed to Appendix A. Note that it completes 
the arguments introduced in Section 10.5.2 of [4], where it is assumed that the F{ are 
continuous. 

It is well known that, in the univariate case, 9 is a constant which does not depend 
on t. In the multivariate case, 9 is required to be a constant on the lines through the 
origin. In the next section, we will take into account this homogeneity property when 
constructing the estimators. More particularly, the multivariate extremal index function 
has the following properties (see [28, 33], Section 10.5.2 in [4], [24] and [25]): 

1. < 0(t) < 1 for all r £ (0, oo) d . 

2. 6(t) is a continuous function of re (0,oo) d and is scale invariant, that is, 9{ct) = 
6(t) for all c> and r G (0, oo) d . 

3. 9 can be extended by continuity to [O,oo) rf \{0}. Let rW = (0, . . . , 0, n, 0, . . . , 0). 
The univariate extremal index, 9i, say, for the ith component sequence {X n ^) n >\ 
exists and 9 i = 9(j( i >). {0} is a discontinuity of 9 if there exist i and j such that 
9i 7^ 9j. Note that the functions H and H can be extended by continuity to [0,oo) d . 
In particular, we have H(t^) = e^, H(t®) = e~ e ' T ' and H(0) = H(0) = 1. 

4. Bounds for 9{t) are given by 

max(6> 1 r 1 , . . .,9 d T d ) < 9^-1 h 9 d r d 

-lnff(r) " -lnff(r) 

The upper bound corresponds to the case where G has independent components and 
the lower bound corresponds to the case where G has totally dependent components. 

5. If G and G have independent components, then 9(t) — Xa=i ^i r i/X)i=i r *- 

In the next sections, we will illustrate our technical conditions and our limiting results 
with three examples of bivariate processes. Let us now introduce these processes and 
discuss their extremal properties. The first process will be considered as the benchmark 
because all the technical conditions can be easily verified and the calculations of the 
asymptotic variances of the estimators can be carried out explicitly. It is the bivariate 
process with independent univariate sequences and independent components. The second 
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process is a bivariate squared ARCH(l) process with independent components. There is 
no cross-sectional dependence, but each component is time dependent. The third process 
is a bivariate autoregressive process of order 1 with dependent innovations. By modifying 
the values of the parameters of this process, we may have cross-sectional dependence or 
independence and temporal dependence or independence. Recall that for d = 2, we have 
(Xi = (Xi tl ,Xi t2 ))i>i, t = (ti,t 2 ) and u n (r) = (u„,i(ri), u„ l2 (r 2 )). 

Example 2.1. The bivariate independent process with independent components: X^i = 
£2.1 and X^ 2 =6,2, where (<!;/, i);>i and (£/. 2 );>i, are two independent sequences of i.i.d. 
standard exponential r.v.s. It is easily seen that u n .i(r) = u n ^{T) = ln(n/r), — ln_ff (t) = 
T\ + t-2 and 9{t) = 1. 

The cluster of exceedances of are of size 1, that is, the cluster-size distribution is 
given by 7rW(l) = 1, ^(k) = for k> 1. 

The associated series is given by zj = max(ri exp(^; j i), t 2 exp(£/ j2 )) and we have 
(x) ~ n>r -1 (n + t 2 ) as n — > 00. 

Example 2.2. A bivariate squared ARCH(l) process with independent compo- 
nents: A ;+M = + A l A / ^)^ 2 +1;i for I > 1 and i = 1,2, where (£i,i)j>i and (£j,i)/>2 
are two independent sequences of i.i.d. standard Gaussian r.v.s, rji > and < 
Aj < 2c Y , where y is Euler's constant. We assume that Xi t \ and Xi t % are drawn 
from the univariate stationary distributions. Let Ki be such that E(Ai^ 2 i ) Ki = 1 for 
i = 1,2. There exist constants Cj such that ^(cc) ~ c;:r~ Ki as x — ► 00. It follows 
that u rh i{T) ~ (nci/r) 1 / Ki as n — * 00 (see, e.g., [21] and [17]). Let = tt{j > 

1 : Xi Y\\ = i{^iQi) > x} where X t is independent of (£j,i)i>i and P(A.; > x) = x~ Ki , 
x > 1, and define p fc)i = P{Ri{l) = k),k>0. 

Since the components are independent, we have — ln_ff (t) = t\ + r 2 and 



firi+f 2 T 2 
0(t) = • 

Tl + T 2 

Moreover, 4 =p ,i = / 1 ° i , (V^iIli= 1 (^) < x)^"^- 1 dx, i = 1,2. 

The clusters of exceedances may be of any size. One can show that the cluster-size 
distribution of is given by 

x = a 77> ni + a 77> n2 ' 

v\T\ + V2T2 OlTi + V2T2 

where Wi(k) = (pk-i,% —Pk,i)/po,i, k>l and i = 1,2 (see [17]). 

Since there is no analytic expression for the stationary univariate distributions, an 
explicit form of the associated sequence cannot be given. 



Example 2.3. A bivariate autoregressive process of order 1 with dependent innovations: 
Xi+i.i = PiXi ti + 6 + i,, for I > 1 and i = 1,2, where (&,i,$,2)j>i is a sequence of i.i.d. 
vectors with a bivariate unit Frechet extreme value distribution, that is, 
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P(Zt,i<xi,&<x 2 )= exp(-( — + — )A[ g )) 
:= exp(-B(xi,x 2 )), 

where A is a convex and differentiable function bounded below by max(i, 1 — x) and 
above by 1. We assume that < pi < 1 and that (Xl.i, Xl j2 ) is drawn from the stationary 
distribution. We have that Fi(x) ~ (1 — pi)~ 1 x~ 1 as x — > oo and it follows that u„.i(r) ~ 
n/((l — /3i)r) as n-> oo. By Theorem 2.1 in [36], we deduce that 

oo 

- luff (r) = ^B(((l - PiKnr 1 , ((1 - pa)^) _1 ). 

fc=0 

Similar arguments as in Section 6 of [33] show that 



e( v = -lnff((ri,r 2 ))+lnff((piri,p 2 r 2 )) 
- lnff ((n,r 2 )) 

It is important to note that if p% = p 2 = p, then — lnff(r) = ^(t^ 1 , t^ 1 ) and 0(t) = 

(1 — p). The multivariate index does not depend on t. If yl = 1, that is, if and £^2 

are independent, then — \\~iH (t) = n + t 2 and 0(r) = ((1 — pxjTi + (1 — /9 2 )t 2 )/(ti + r 2 ). 

(V) 

The clusters of cxcccdances may be of any size. The asymptotic distribution of iV„ 
may be obtained by using results of Section 2 in [6]. Moreover, since there is no analytic 
expression for the stationary bivariate distribution, an explicit form of the associated 
sequence cannot be given. 



3. Defining the estimators 

In this section, we explain our approach to estimating the multivariate extremal index 
function. As in [37], we consider a block declustering scheme and estimate intermediate 
thresholds such that we only have to take into account the block length to study the 
asymptotic distribution of the estimators. 

Let us divide [1, ...,n] into k n blocks of length r n (k n is the integer part of n/r n ), 
Ij = [(.? - l)r n + 1, . . . , jr n ] for j = 1, . . . , k„ and a last block h„+i = [r n k n + 1, . . . , n] . 
The number of exceedances for the jth block is defined by N^- = J2ieij ^{x t ^u r (x)} 
for j = 1, . . . , k n , where it rnj j(r,) = F?~~ (1 — Ti/r n ), i = 1, . . . , d. The main issue when using 
these quantities to construct estimators is that the thresholds u rn ^(ri) are unknown since 
they depend on the univariate marginals of the stationary distribution. They have to be 
estimated from the data. As in [37] , we consider estimators of the thresholds which are 
based on the order statistics. If < r < r n , let u rniJ (r) = X^ knT ^ i, where -XVfcVi is the 
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fcth largest of X^j, . . . , Xk n r n ,i and \x\ denotes the smallest integer greater than or equal 
to x. If t = 0, let u r „,i(0) = oo. Now, define N£'j = l{Xi^a rn (x)} f° r T e [0,?"n] d . 

In order to estimate the multivariate extremal index function, it seems natural to 
exploit the characterization given by (1.4). Let = iVn ((0, 1]) and note that, under 
appropriate conditions (see the following section), 



lim P(N^ = 0) = lim P(M n < uJt)) = H(t) 

n — >oo n — >oo 

and that, by (1.2), 



lim E(ArM)= lim n(l - F(u n (r))) = -InHM. 

n — >oc n — >oc 

Let us use the empirical distribution of the number of exceedances to provide empirical 
counterparts of H(t) and — IiiH(t). We define 



^™( T ) = I^EW> =0} and -^H n (T) = ^f]N% j forre[0,r„] d . 

One may consider — \nH n {r)J{— hiH n (T)) in order to estimate 9(t). But, unlike the 
multivariate extremal function, this function is not scale invariant (see Property 2 in the 
previous section). Hence, we introduce a first estimator which satisfies the homogeneity 
property: 



C)(r) = -lnfln(r/I(r)) | ^ g 
-\nH n {T / L{t)) 

where L is a known function from [0,oo) d \{0} to (0,oo) which is homogeneous of order 
1. For example, consider the family L C:a (r) = c(^\ =1 Ir^ ) 1 / for a > and c > 0. 

The second estimator is derived from the characterization of Proposition 2.1. Let us 
consider the number of exceedances of (Z^) n >i above the threshold v n T \x): 



n 

i=i 



Proposition 2.1 implies that lim n ^oo — lnP(iV„ = 0) = 9{t)h. In order to construct 
an alternative estimator of the extremal index function, we can follow the approach 
developed in [37]. First, we replace the z\ by their empirical counterparts since the 
marginal distribution functions Fi are unknown. Let Ri t i denote the rank of Xu among 
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(X\ t i, . . . ,Xk n r n .i)- In the case of ties, the lowest rank for the ties is used for each tie. We 
define 



where 

k n T n 



Li 



kn^ri ~t" 1 Rl,i 



We then introduce the number of exceedances of {Z\ for the jth block: 

Eig/- ^{z (t) >^ t) (x)}' ^ s P rev i° us ly; w» ( x ) i s unknown. However, it may be esti- 
mated by vi^\>c) = ^-ly where Z^ is the ([/s„>r~])th- largest value among 
Z^ , . . . , ^ 7 "^, . Finally, let us define N^'J' as the counterpart of N^'J' , where Vr„ (x) 
is replaced by ir^\>c), and introduce the second estimator 



9® (r) = -x- 1 In ^i- £ l {i v( :: j) =0} ^ , t G [0, oc) d \{0}, x £ (0, r n ]. 
Note that this estimator is scale invariant without transformation on t. 



Remark 3.1. In [37], three estimators of the univariate extremal index are introduced. 
The first estimator, denoted by §[ \ , is very close of our estimators 0^ and when 
they are evaluated at the points r = t", i = 1, . . . , d. In fact, if L is assumed to be a 
constant equal to 1 and >t = tj , we have 



^(^"©p^fri and W'MS?- 

It follows that in the univariate case (i.e., d = 1), both estimators have the same asymp- 
totic behavior as 



4. Main result 

In this section, wc first present and discuss technical conditions which arc required for the 
asymptotic normality of the estimators. These conditions arc quite similar to conditions 
introduced in [37] which are used, in particular, to establish the asymptotic properties of 
the estimator of the univariate extremal index §[ T ^ (see Remark 3.1 above). They might 
appear quite stringent in comparison with those of [45] , where the asymptotic properties 
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of the blocks and runs estimators of the univariate extremal index are studied. This 
is not the case for two reasons. First, we estimate intermediate thresholds and do not 
consider them as tuning parameters, contrary to [45]. This allows us to establish the 
asymptotic properties of intermediate empirical processes, which is more complicated 
and necessitates more conditions. Second, these conditions guarantee the existence of the 
asymptotic variances of the estimators, whereas [45] just assumes the convergence of the 
variance of a partial sum to the asymptotic variance and does not give any condition such 
that this convergence holds. Finally, one can refer to Section 4 of [37] for a comparison 
of similar conditions to those in [38] that are needed for convergence of the tail empirical 
process of a univariate stationary sequence. 

Let us turn to some definitions which are the natural multivariate versions of definitions 
from [32] (see also [28] and [33]). 

Definition J^.l. Fix an integer m > 1. Let T p%q = -Fp, g (Ti, . . . ,T m ) be the a-algebra 
generated by the events {Xi < u n (Tj)}, p <l < q and l<j<m, and let 

a n ,i{n, r m ) = sup|P(A (IB) — P{A)P(B) :AeF u ,Be F t+l , n , \<t<n-l\. 

The A({u ra (Tj)}i<j< m )-condition is said to hold if limn^oo a n j„ (ti, . . . , r m ) = for 
some sequence l n = o(n). 

Definition Suppose that the A({u n (Tj)}i<j< m )- condition holds. A sequence of 

positive integers (q n )n>i is said to be A({u„(rj)}i<j< m )-scparating if, as n—> co, q n — 
o(n) and there exists a sequence (l n )n>i such that linin^oo nq~ x a n _\ n (ti, . . . , i~ rn ) = 
and l n = o(q„ ). 

We now give a decomposition of the numbers of exceedances when considering two 
vectors of thresholds, u„(ti) and u n {T2) for ti,t 2 G [0,oo) d . We define 

p 

/vr(Tl,-r 2 ) _ \ ~* -, 

iV n,0,p — 1 {X i ^u„(-ri)}U{X,^u ?l (-r2)}' 
1=1 

P 

A r(-r 1 ,-r 2 ) _ \ " -, 

1Sl n,l,p ~ i {X,^„(T 1 )}\{X 1 ^„(T 2 )}. 
(=1 

P 

at(ti,t 2 ) _ -, 

JV ",2,p - 1 {^^n(T 2 ))\{X,^„(T 1 )). 

(=1 

p 

n (ti,t 2 ) 

ly n,3,p - 2^ i {X,^ M „(-r 1 )}n{X ! ^«„(-r 2 )}- 
(=1 

Note that N { n 2' p T2) = E?=i and = + N^\ i = 1,2. 

We continue by presenting the first technical condition and then discussing the weak 
convergence of the sequence (N^'^ , Nfy£ 2 \ ) n >i- 
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Condition (CI). 

(i) The stationary sequence (X n ) n >i has a multivariate extremal index function 8 > 

0. 

(ii) For each t\,t<x G [0, oo) d \{0} ; the A(m„(ti), tt„(r 2)) -condition holds and there 
exists a probability measure jt^ 1 ^ 2 ) such that for all ii >0,i2 > 0, 13 > 0,ii + 12 -M3 > 1, 

tt ( ^ 2) * 2 , k) = n Um ^(^;;„ 2) = 4, ft = 1, 2, 3|< V 9 ; 2) > 0) (Cl.a) 

/or some A(m„(ti), n. n (r 2)) -separating sequence (q n )n>i- 

For Tj = (tis, . . .,T d ,i) and i = 1,2, let n V r 2 = (n,i V ri >2 , • . .,r d) i V r dj2 ). For each 
r 1) T 2 G [0,oo) d \{0}, let C = (& 1 1 ' T ' ) , Cfi 1 '^ 5 . Cj,3 1,Ta) )l>i be an sequence of i.i.d. vec- 
tors of integer r.v.s with distribution n^ 1 ^ 2 ^ and 77(1-1,1-2) be an r.v. with Poisson 
distribution and parameter —0{t\ V T2) ln(ff (n V r 2 )) independent of the sequence £. 

The probability measure tt^ Ti ' T2 '> is a key parameter to characterize the distribution of 
the limiting two- level exceedance point process (see Theorem 2.5 and its proof in [32] for 
the univariate case). The following proposition is concerned with the weak convergence 
of the related sequence of the numbers of exceedances. 

Proposition 4.1. Suppose that (CI) holds. Then, 

/At(ti,T2) at(ti,T2) »t(ti,T S )\ D , ^(ti ,T 2 ) W (ti,T 2 ) jir(Tl,T 3 )\ 

^ iV n,l,n ' Jv n,2,n ' Jv n,3,n / ^l Jv l ) Jv 2 ' Jv 3 J , 

(4.1) 

l)(Tl,T 2 ) 

O \ " ^(ti,t 2 ) >(ti,t 2 ) >(ti,t 2 )\ 
— ''i,2 >M,3 /• 

2=1 

Moreover t^ 1 " 1 ' 1 " 2 ) is scofe invariant, that is, for each Ti,t 2 G [O,oo) d \{0} and c > 0, 

^(c-riXT-2) _ ^(1-1,1-2) _ 

The proof of Proposition 4.1 is postponed to Appendix A. 

Let us examine the distribution of the cluster sizes of the two- level exceedances, 
7r (Ti,T 2 ) ) f or tjjg £ rg ^ exam pi e introduced in Section 2. The distribution for the second 
example may be derived by considering its Laplace transform. The distribution for the 
third example may be derived by using results of Section 2 in [6] . 

Example 2.1 (continued). The clusters of the two- level exceedances are of size or 
1. More precisely, the distribution of the cluster sizes is given by 

^•^(1,0,0) = (^■i-^.2) + + (ri,i-r 1 , 2 )+ ) 

Tl,l V 71,2 + T2.1 V T2, 2 

rr^^)(0,l,0)= (r2 - 2 ' r2 - l)+ + (Tl ' 2 ' ri - l)+ , 
n,i Vri,2+T2,i Vr 2 , 2 



1038 C. Y. Robert 

^ (Ti ' T2) (o,o, i) = i - 4 T1,Ta) (1,0,0) - 4 T1,Ta) (o, 1,0). 

Let us turn to the second technical condition which is a multivariate version of Con- 
dition (C2) in [37]. Note that, since the estimated thresholds for our estimators are con- 
tingent on k n Ti, i = 1, . . . , d, and k n *c, and since k n may be chosen up to a proportional 
factor, we can assume, without loss of generality, that r and k are bounded. Hence, let 
us now assume that r G [0, l] d . 

Condition (C2). 

(i) Let r > 4d. There exists a constant D = D(r) > such that for all T\, T 2 G [0, l] d , 

supE^^ -N£')\ r <D\n -r 2 |. (C2.a) 

(ii) Let lu > (Ad — l)r/ (r — 4c?). There exists a constant C > such that for every 
choice of Ti, . . . , T m G [0, l] d , m > 1 and n > I > 1. 

a n ,i(ri,...,T ro )<a;:=Cr w . (C2.b) 

(hi) (r n ) n >i is a sequence such that r n — > oo anc? r„ = o(n) and there exists a sequence 
(l n )n>i satisfying 

l n =o{rl /r ) and lim nr^a^ = 0. (C2.c) 

n — KX3 

Let us describe some intuitions regarding this condition. First, (C2)(i) restricts the 
size of clusters by assuming that Nn ~ has a suitably bounded rth moment. It 

provides an inequality which will very useful to prove tightness criteria for intermediate 
empirical processes introduced in Section 5. Note that 



i=l 

It follows that it is sufficient to show that for each i = 1, . . . , d, there exists a constant 
A > such that for all n,T 2 G [0, l] d , 

S upE(4 T?>AT = ,)) -ivf^'Y < Aki-7?, 2 |. (4.2) 
n>l 

It is easily seen that (C2)(ii) is satisfied by strongly mixing stationary sequences where 
the mixing coefficients vanish with at least a sufficient hyperbolic rate. The underlying 
idea of the block declustering scheme is to split the block Ij into a small block of length 
l n and a large block of length r n — l n . (C2)(ii) and (C2)(iii) essentially means that l n is 
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sufficiently large such that blocks that are not adjacent are asymptotically independent, 
but does not grow too fast so that the contributions of the small blocks is negligible. Let 
us give now some clues explaining why this condition holds for the examples introduced 
in Section 2. 

Example 2.1 (continued) . {{X^i,Xi^))i>i isani.i.d. sequence. Therefore, a n ,z(i"i, • ■ ■ , 
r m ) = for every choice of Ti,...,T m G [0, l] d , m > 1, n > 1, I > 1. Moreover, 

( x W Ax W) ( T (0 VT {i)) 

iV n 1 2 — Nn 1 2 has a binomial distribution with parameters n and |r,.i — T^l/ft- 
Condition (4.2) is easily verified for any integer r. 

Example 2.2 (continued) . The components of {{Xi t i,Xi^))i>i are independent and 
each component is geometrically strong- mixing (see Example 3.1 in [37]). Moreover, 
bounds for the moment condition (4.2) can be obtained for any integer r by using the 
same arguments as for Lemma 6.1 in [37]. 

Example 2.3 (continued). {(Xi t i,Xi } 2))i>i is a bivariate positive Harris recurrent 
Markov chain. Moreover, it is a particular case of a first-order stochastic equations with 
random coefficients. Following [26], one can show that {{Xi^i,Xi^))i>i is geometrically 
absolute regular and strong-mixing. Moreover, bounds for the moment condition (4.2) 
can be obtained for any integer r by using the Markov property of the components and 
similar arguments as for Lemma 6.1 in [37]. 

Finally, we assume that the convergence rate of r„ to infinity is such that the bias 
of our estimators is asymptotically negligible with respect to their variance. Moreover, 
we need a condition on the regularity of H and H to guarantee that the asymptotic 
distribution is Gaussian. 

Condition (C3). 

(i) The sequence (r n ) n >i satisfies 




and 

lim yfk~ n sup \r n (l-F(u r (r))) + InH (r)| = 0. 

"- 00 xe[o,i] d 

(ii) The functions H and H are (Frechet) differ entiable on (0, l) d and their deriva- 
tives can be extended by continuity to [0, l] d . 



Example 2.1 (continued). Note that as n 



^ sup \P(N { r l\=Q)-H(T)\~ e - 
-re[o,i] d " 



3-i nV2 

2 3/2 : 



1040 



C. Y. Robert 



k n sup \r n (l-F(u rn (T)))+lnH(r)\ = ^— 
-re[04] d r n 



It follows that if r n =o(ri 1 / 3 ), then Condition (C3) does not hold. 



We end this section by giving the distributional asymptotics of the estimators. Let 
O(-) be a pathwise continuous Gaussian process on [0, l] d \{0} with covariance function 
given in Appendix B. 

Let V l = {t:t/L(t) g [0,l] d \{0}} and 4^ = {t : xt/(- lniJ(r)) e [O,l] d \{0}}. 

Theorem 4.1. Suppose that (CI), (C2) and (C3) hold. If we letm>l and n, . . . , r m £ 

n/^C(tO - e{r t )) i=l _ m Z (e^L- 1 ^)))^!,...,™. 

// we let m > 1 tmrf Ti, . . . , T TO G i/ien 



vC(^ 2) (n)-^(T i )) j=li ... )m -(e(T i x(-ln^(T,))- 1 )) i=li ... im . 

Although the estimators are very different from the point of view of their construction, 
they share the same asymptotic distribution up to a proportional factor. 

Example 2.1 (continued). The calculation of the asymptotic variance of §£\t) and 
0u\t) can be carried out explicitly. Let M(r) = — In i7(r)L _1 (r). We have 

Var(e(xL- 1 (x))) = M(r)- 2 (e M W - l - M{r)), 
Var(6(r^(-ln#(r))- 1 )) = (x^e 3 * - 1 - x)). 

It is worth mentioning that the asymptotic variance of 6\ (t) does not depend on r. It 
is smaller than the asymptotic variance of e { n\r) if x <M(r). 

Note that if r = r' 1 ' , _L = 1 and x = tj, then we obtain the same asymptotic variance 
as for 8i~n ( see Remark 3.1). 

Example 2.2 (continued). Let us characterize the asymptotic variance of the second 
estimator. We have 

Var(e(r^(-lni7(r))- 1 )) 

-J ( (hn + e 2 r 2 \ e in + e 2 r 2 \ 

= x exp x — 2x 1 

V V T 1 +T 2 J Tl+T 2 J 
+ X 

n + T2 
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(717 1 1- (7 2 7 2 S ' I l*\ ( * I 1\ ( (T (1> T {2 >) , (x (2) T (1, )V • 7 \ 

l6>2 ; > (« + «)C? +0( 7lA ; +7T^ ' x ^)(l,J,fc). 

>_0,j>0,k>< 
i+j+k>l 



Ti + To 

i>0j">0,fc>0 



For the first estimator, replace >e by -Ajf(r). A comparison between the two asymptotic 
variances is not obvious. 

Note that if r = , L = 1 and x = T il then we obtain the same asymptotic variance 
as for (see Remark 3.1). 

It is possible to weaken Condition (C3)(ii) by assuming that there exists an open set 
O included in (0, l) d where the functions H and H are (Frechct) differcntiable. One then 
has to replace * L by {r : t/L(t) G O n [0, l] d \{0}} and by {r : xt/(- lniJ(r)) e 
On [0,l] d \{0}} in Theorem 4.1. 



5. Simulation study 

In this section, a simulation study is conducted to investigate the performance of the 
estimators on samples of moderate size. Data are simulated from: 

• the bivariate independent process with independent components of Example 2.1 - 
we have 6(t) = 1; 

• the bivariate squared ARCH(l) process with independent components of Example 
2.2 - we choose ry = 2 x 1(T 5 , Ai = 0.7, A 2 = 0.3 so then we have (sec [17]) 



„. , 0.579ri + 0.887r 2 
= ; 

• the bivariate autoregressive process of order 1 with dependent innovations of Exam- 
ple 2.3. We choose p\ = p 2 = 1/2 and B(x\,X2) = (%i 2 + x 2 2 ) 1 ^ 2 so then we have 
0(t) = 1/2. 

We study the performances of the first estimator 6^ associated with the functions 
Lc.aij) = c(X)j = i l''";!' 1 ) 1 ^' 1 for c = 2 and a = 1, c = 1 and a = 1, c = 2 and a = 2, c = 1 and 
a = 2, and we undertake comparisons with the second estimator On associated with x = 
1. For each process, we generate 500 sequences of length n = 2000 and for each sequence, 
we compute the estimates for r = (cos </>, sin cf>) with </> = kn/22 and k = 1, . . . , 10. 

Figures 1, 2 and 3 show the means (left) and the root mean squared errors (RMSE) 
(right) of the estimates as functions of the angle, <fi, for k n = 50, 100, 150, 200 and for the 
three processes. First, observe that the bias of the estimators decreases as the size of the 
blocks increases. The estimators are nearly unbiased for k n = 50 and k n = 100, but they 
show a positive bias when k n = 200, except in the case of the bivariate squared ARCH 
process and for the large values of (f>. Conversely, the variances of the estimators increase 
as the size of the blocks increases. This is the ordinary variance-bias trade-off encountered 
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Figure 1. The i.i.d. sequence. Left: means of the estimates of the multivariate 
extremal index function; the gray solid line represents the true function. Right: 
RMSE of the estimates of the multivariate extremal index function. The estima- 
tors which are considered are §„ associated with Z/2,1 (- -), (• ■ • ■), 
Ij2,2 (- - - -), £1,2 ( ) and 8^ associated with >c=l ( ). The graphs show the aver- 
age over 500 samples. 



with the blocks estimators. Note that the minimum of the RMSE is generally observed 
for large values of k n (150 or 200). 
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Figure 2. The bivariate squared ARCH(l) process. Left: means of the esti- 
mates of the multivariate extremal index function; the gray solid line repre- 
sents the true function; Right: RMSE of the estimates of the multivariate ex- 
tremal index function; the estimators which are considered are 0„ associated with 
I/2,i (- — ), Li,x (■ ■ ■ ■), 1/2,2 (- - - -), Lx,3 (- - -) and §„ associated with h = 1 
( ). The graphs show the average over 500 samples. 



For the i.i.d. sequence, the first estimator associated with the function £2,1 performs 
uniformly better. The reasons for this may be that the asymptotic variance is smaller 
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Figure 3. The bivariate autoregressive process. Left: means of the estimates of the multivariate 
extremal index function; the gray solid line represents the true function. Right: RMSE of the 
estimates of the multivariate extremal index function. The estimators which are considered are 

8^ associated with Z/2,1 ( ), Li,i (• ■ ■ •), £2,2 (- - - — ), £1,2 ( ) and 8^ associated with 

x=l ( ). The graphs show the average over 500 samples. 



and that the estimated thresholds are higher than those with the choice c = 1 and hence 
the bias is smaller. The RMSE is mimimal for k n = 200. 
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Figure 4. The i.i.d. sequence. Ratios between the sample variances and the asymptotic vari- 
ances for the estimators of the multivariate extremal index function are shown. The estimators 

which are considered are 8^ associated with Z/2,1 ( ), L\,i (• • • •), 1/2,2 (- - - — ), £1,2 (- - 

-) and ^ associated with x=l ( ). The graphs show the average over 1000 samples. 



For the bivariate squared ARCH(l) process. On associated with the function L\ 1 
and 8„ perform better than the other estimators. Note that, as for the previous process, 
— InH and L^i are equal, which explains why both estimators have the same asymptotic 
variance. The RMSE is mimimal for k n = 150, except for the large values of <fi. 

For the autoregressive process, there is a relatively small sensivity of the estimates to 
the choice of the estimator when k n = 200. The second estimator always performs better. 

Overall, 6 n appears as a good candidate to estimate the multivariate extremal index 
function. Its performance on samples of moderate size is often better than the perfor- 
mance of 9n^ and, moreover, it does not necessitate the choice of a tuning function. 

Figure 4 shows, for the first process, the ratios between the sample variances and the 
asymptotic variances for the estimators of the multivariate extremal index function. It 
illustrates that for sequences of length at least n = 2000, the variances of the estimators 
can be well approximated by the asymptotic variances when they can be calculated or 
estimated. 



6. Intermediate results and proof of Theorem 4.1 

The proof of Theorem 4.1 and some results related to the weak convergence of interme- 
diate empirical processes are gathered in this section. We let K be a generic constant 
whose value may change from appearance to appearance. 
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We first introduce the Skorokhod space of ladcag multiparameter functions and give 
the essential ingredients that will be required for characterizing the asymptotic behavior 
of the intermediate processes. 

Let Bd be a cube in M. d . If r £ Bd and if. for i = 1, . . ., d, Ri is one of the relations < 
and >, then let Qi? ll ...,/j d (T) be the quadrant 

{a- = (oi, . . . , aa) £ B d : OiR^i = l,...,d}. 

We denote by D(Bd) the space of functions from Bd to E which are "continuous from 
below, with limits from above" in the sense defined by [1]. More precisely, / £ D(Bd) 
if, for each r £ Bd, Jq(t) = lim CT _>T- j(Te Q /(er) exists for each of the 2 d quadrants Q = 
Qi?i,....i?. d ( T ) and /(t) = Jq < <. Let us assume that it is equipped with the metric d° 
which is equivalent to the Skorohod metric d and such that it makes D(Bd) a complete 
separable metric space (see [1], Section 2 and [3], Section 12). 

A sequence (f n )n>i of D(Brf)-valued processes converges weakly in the Skorohod topol- 
ogy to a £'(B c i)-valued process / (/«(•) => /(•)) if ^f(fn) — > E</2(/„) for all Skorohod- 
continuous bounded functions ip:D(Bd) — > R. A criterion for the weak convergence of 
D(Bd)-valucd processes can be given in terms of the weak convergence of the correspond- 
ing finite-dimensional distributions together with a tightness condition (see Theorem 1 
in [1] and the proof of Theorem 6.1). 

It is often convenient to consider the restrictions of the functions of D{Bd) to a subset 
of Bd- If Cd is a cube included in Bd and if / £ D(Bd), we denote by rc d f the restriction 
of / to Cd- We have the following convergence property: if /„(•) => /(•) in D(Bd) and / 
is continuous at the lower boundary of Cd, then rc d f n (') fc d f(') m D{Cd) (see, e.g., 
Lemma 4.17 in [35] for the univariate case). 

We now turn to the definition of the intermediate processes and characterize their 
asymptotic distribution. First, let us introduce, for r£ [0, l] d , 

V n (r) = Vh t (H n (T) - P(N£> A = 0)), 
W n {r) = \fK{{~ lnff„(r)) - r n P(X l £ u rn (r))), 

where 

ff «(T) = ^-Ev- ) .=o} and - lni w = f E< x) r 

3=1 3=1 

We define the first intermediate D([0, l] d ) x D([0, l] d )-valued process by U n {r) = 
Wn(r))'. Observe that U n depends on the unknown vector of thresholds M rn (x) 
and cannot be used in practice. In the univariate case, W n is called the tail empirical 
process and has been studied for dependent sequences in [7, 8] and [38]. 

Theorem 6.1. Suppose that (CI) and (C2) hold. There exists a pathwise continuous 
centered Gaussian process U with covariance matrix C(-,-) = {Ci,j{','))i<i,j<2 given in 
Appendix B such that U n (-)=>U(-) = (V(-), W{-))' in D([0,l] d ) x D([0,l] d ). 
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The proof of Theorem 6.1 is presented as a series of two lemmas. Let us define the 
large blocks Ij and the small blocks /* by, for j = 1, . . . , k n , 

Ij = [(j - l)r n + 1, • ■ -,jr n - l n ], I* = [jr n - l n + 1, . . .,jr n ]. 
We introduce the quantities 

ff n ( T ) = r L £ 1 {Jv w -. A =o}' ff »W :; T^ 1 c w ' 4 ^ |Tl, *>»}' 
-ln^(r) = ^E^lf, -lnff:(r) = ^IX^*, 

3=1 J = 1 

and consider the following processes 

^„ A (r) = ^{H£{r) P(N^i A = 0)); 
K(r) = VK{K{t) + P(N^f = 0, N%>f > 0)); 
W*(t) = vC((-lni/,f (r)) - (r„ - Z„)P(X ; ^ ^„(t))); 
W„*(t) = v^((-ln^(r)) - Z„P(X ( £ u r „(r))); 

= (^ A (r), W*(t))', C/;(r) = (y*(r), W* n {r))' . 

Note that H n (r) = H^{t) + H*(r), -\nH n (r) = -\hH£(t) - ln#*(r), V n (r) = 
V n A (T) + V:(t), W n (r) = W r ?(r) + W*{t) and U n (r) = C/ A (r) + LT(r). 

Lemma 6.1. Suppose that (CI) and (C2) hold. Let m>l and ti,...,t to € [0, l] d . 

Then, 

(C7 n (r i))i=i ) ..., m —> {U{t i))i=i,... )TO . 
Proof. By similar arguments as in Lemma 6.6 of [37] with 

2lj 

Ad V < v < r 

lu-1 

(this can always be assumed), we have U n (r) 0. It is only needed to be checked that 

(t7^(T i )) i= i,..., TO ^(E7(r i ))<=i,...,m. 
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But we can use the arguments of Lemma 6.7 in [37] and replace (C0)(b) in [37] by 
(C2)(a) with t 2 = in order to establish the weak convergence and to conclude that 
(U(Ti))i—i t ,,, !Tn is a Gaussian centered random vector with covariance matrix 



Gov(U(T l ),U(r k ))=E(U{Ti)U'(T k ))=C{r l ,T k ), l<l,k<m. □ 

Lemma 6.2. Suppose that (CI) and ( G2) hold. Let us define the modulus of continuity 
of f e D([0,l} d ) by 



w f (S) = su P {|/(x) - f(r')\ :t,t' e [0, l] d , |r - r'| < 5}. 
Let e > 0. Then, 

lim lim sup P(w Vn (S) > s) = 0, (6.1) 

hmlimsupP(unv„(<5) > e) = 0. (6.2) 

Proof. We combine some arguments from Section 5 of [29] and some arguments from 
the proof of Theorem 1 in [5]. Let L(2 TO «) be the set of all points (h,.. . 1 l d )/2 m " with 
[ i G{0 J l,...,2 m -} I i=l,...,d. 

Since H n {r) is a monotonically non-increasing function in each component of r and 
— In H n (r) is a monotonically non-decreasing function in each component of t, we have 
(see [31], page 262) 

w v JS) < 6su P |K(ti) - V n (r 2 )\ +4dy*~2- m ", 
wwM < esupl^^xx) - W„(t 2 )\ +4d^2- m ", 

where the "sup" is to be taken over all ti,t 2 E L(2 m ") with |n - r 2 | < 5 + 2~ mn+1 . If 
m„ is chosen such that linin^oo \fh^2r ra ' n = 0, (6.1) and (6.2) will follow if we can show 
that 

limlimsupP(sup|K(Ti) - K(r 2 )| > e) = 0, 

8->0 n 

limlimsupP(sup|W n (Ti) - W n (r 2 )\ > e) = 0, 
where the "sup" is to be taken over all Ti,t 2 G L(2 mn ) with |xi — t 2 \ <S. Let us define 

Y Vli (T lt T 2 ) = (1^=0} - P ( N r T J = 0)) - (1^3=0} - P (^S = °))> 

Y Wj (t u t 2 ) = (N^j - r n P(X l £ Ut b (ti))) - (Jv£J - r n P(Xi £ u rn (r 2 ))), 
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fe„ fc„ 
<SV,„(ti,t 2 ) = y^Y v>j (T 1 ,T 2 ), S W>n (T 1 ,T 2 ) = ^2y wj (t 1 ,t 2 ). 

3 = 1 3=1 

We now want to use equation (4.3) of Theorem 4.1 in [39]. Note that a rrl i < cti for 
I > 1 and r n > 1. Let 2<v<p<r< oo, n > 0, and assume that u> > v/(v — 2) and 
u > (p — l)r/(r — p). We choose p — Ad and r = 2di> with 2 < w < 4, which leads to 
uj > (Ad - l)r/(r - Ad). We deduce that 

E|5v, n (r 1 ,T 2 )| 4d <ei|iV J -(ri,T 2 )||^ + k 1 +»\\Y V j(T 1 ,T 2 )\\% v , 
ElSwinCri.raJl^^^liy^CTx.raJII^ + ^liy^frx.ra)!!*!,. 

Note that, for A > 1, 

|y^-(Ti,T 2 )i A < 2 A (|i {JV ( :: , =0} - i {JV (-) = o } i A + l p (<5 = °) - p ( N r?J = °)i A )' 

|ywu(Ti,r a )| A < 2\\Nl T J - N^\ X + \EN%$ - EN { r T J \ x ). 
Since 

we deduce by (C2)(a) that for 1 < A < r, 

E l 1 {<^=o} - V; 2 >=o} |A - m ^ ~ N ^ x - DlTl T2 '- 

Moreover, 

\P(N£J = 0) - P(N™ = 0)| < n\ N ^ =0} l { ^)=o } l < E l<^ - N r T J\ 

and 

|EiV r y - EivQ | < E\N^J - I < E\N™ - iV^ \ r . 

It follows that for A > 1 and \t\ — t 2 \ < 1, 

flEJV^ - E7V r ( ^|) A < (E|iV^V - <^D A < (Din - r 2 |) A < if|n - r 2 |. 

We deduce that for |ti — T 2 | < 1, 

\\YvArur 2 )C < K\ Ti -t 2 \^, \\Y V j{r lt T 2 )\\^ < K\r x - t 2 \ 2 '\ 
\\Y w AruT 2 )\t d < K\n - t 2 \^/\ \\Y w , 3 {T U T 2 )\\f dv < K\n t 2 \^ 
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and it follows that for any n > 0, 

E\S v , n (T U T 2 )\ M < K(k 2 n d \ Tl - r 2 | 4 ^ + fci+ K |n - t 2 | 2 /"), 

E|S w ,„(n,T 2 )| 4d < x(/c^|ri - r 2 \ 4d ^ + fci+«|Ti - r 2 | 2 /' u ). 

Since right-hand sides of the previous two inequalities are the same, we only consider the 
case of the process W n . We then deduce that 

E|W„(ti) - W n (r 2 )\ M < KQn - r 2 \ Ad /» + kl+ K - 2d \ Tl - t 2 \ 2 ' v ). 

If fci+ K " 2d < |n - T2 \ id l' v ' 2 l v or, equivalcntly |n - t 2 | > tf^+x-WtWd-D) ^ thcn 

E|W„(n) - ^„(r 2 )| 4d < K\n - T 2 \ id '\ 

In particular, if t 2 =n + e,;2~ 7 , where e,; = lW such that 2 ' > fc„ +re 2d )/( 2 ( d_1 ) ; we 
get 

E|W„(n) - W„(r 2 )| 4d < K(2- r ) 4d A'. 

Let m(<5) = max{7 € N:52 7 < 1} and < a < 1. By using the same arguments as in 
Section 5 of [29], we have 

P(sup{|W^(n) - W n (r 2 )\ :ti,t 2 G L(2 m "), |n - r 2 | < 5} > e) 

d m n 2 7 2^—1 2 7 

<E E E---E---E p (i^y)-^y + ^ 2 ~ 7 )i 

i=l 7=m(<5)3i=l ji=0 j d =l 

>(l-a)a~<- m ^E(4 : d 2 )- 1 ) 

where j = (ji,...,j d ) G i(2 7 ). If 2""? > 2~ m " > j£(i+*-a«0/(2(2cl-i)) j then we ^ by 
Chebyshev's inequality, 

P(\W n (j) - W n (j + e t 2^)\ > (1 - a)a 7 -" i W £ (4d 2 )- 1 ) 
< E|W n (n)- W n (r 2 )| 4d < tf(2-T) 4d /" 



((l-a)a''- m ( 5 )e(4d 2 )- 1 ) 4d " ((1 - a)a'r- m ( s h{4d 2 )- 1 ) id ' 
It follows that 

P(sup{|W n (n) - W n (r 2 )| :n,r 2 G L«(2 m "), |n - r 2 | < 5} > e) 

< ^ V (o-1\d{i/v-V, 1 

" ((1 - a)e(4d 2 )- 1 )4 d j^ 4) a7-"»W 
= d_K ^ 1 

■7 — 771(5) 
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Let us choose a such that 2 d ( 4 /" _1 )a > 1, which is possible since v < 4, and let us choose 
m n such that 

lim v/A~2" m " = and Urn jt«(i+"-a«0/(2(ad-i)) 2 m B = Q) 

n — >oo n — »oc 

that is, such that 

fc -<i/2-«/(2(2d-i))) =o(2 - m „ ) and 2 -" l "=o(fc I ; 1 / 2 ). 

This is clearly possible since v > 2 and k is arbitrarily small. Now let m„ tend to infinity. 
The infinite series converges since 

2 d(4/v-i) a > 1 Fi na ii y) i et § ten( j to q orj C q U ivalently, 
m(i5) tend to infinity. The upper bound tends to zero since 

2 d(4/v-l) > Y and the regult 

follows. □ 



Proof of Theorem 6.1. Lets define the diameter of a rectangle as the length of its 
shortest side. Call a partition of [0, l] d formed by finitely many hyperplanes parallel to 
the coordinate axes a S-grid if each element of the partition is a "right-closed, left-open" 
rectangle of diameter at least S and define w'^ ^(5) : D([0, l] d ) — > R by 

w'f (S) = inf max sup |/(r) — /(<r)|, 

where the infinimum extends over all 5-grids A on [0, l] d . Let us define lis : D([0, l] d ) — > 
R 5 by Il5(/) = (f(s)) seS for each finite set S C [0, l] d . Let T be the collection of subsets 
of [0, l] d of the form Ui X • • • X Ud where each Uj contains and 1 and has countable 
complement. According to Theorem 2 in [1], V n ^>V (resp., W n W) if and only if 

(i) Us (14) Hs(V) for all finite subsets S of some member of T (resp., ns(W„) 

n s (W0); 

(ii) hm^-fO limsup n P(w' v (S) > e) = for all e > (resp., lim^o limsup^ P{w' v (S) > 
s) = Q). 

By Lemma 6.1, we derive the first condition. Now, according to equation (1.7) in [29], 
we have 

w' Vn (S)<w Vn (2S) and w' Wn {8) < w Wn (26), < 6 < 1/2. 

By Lemma 6.2, we derive the second condition. Moreover, by using the same arguments 
as in the proof of Theorem 15.5 in [2], we can show that V and W belong to C([0, l] d ), 
the subset of D([0, l] d ) consisting of continuous functions. □ 

We now substitute the unknown vector of thresholds in the first intermediate process by 
its estimate, u rn (T), and replace P(N^ j = 0) and P(Xi ^ u rn (r)) by their respective 
limits. Let us introduce 

V n (T)=y/K(H n (T)-H( T )), 



W n (r) = y/K((- \nH n (r)) - (- lnff(r))) 
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and define the second intermediate D([0, l] d ) x D([0, l] d )-valued process by tJ n {r) 
(V n (r), W n (T))' . We now establish the weak convergence of this process. 

Proposition 6.1. Suppose that (CI), (C2) and (C3) hold. Then, 
in D([0, l] d ) x D([0, l] d ), where 

T-Tt \ T T / \ ( -VFM'ZfT) 

^ T ) = ^ + U-Hr)Vff(r)'Z(r) 

Vtf (r) = (dH(T)/dT t ) l=1 _ d , V#(t) = {dH{T)/dri) i=lt ... 4 and Z(t) = (VF(^(r))) l=1 ,... id 
with 7Tj(r) = t™, i = 1, . . .,d. 

Note that [7 is well defined on [0, l] d \{0} and can be extended by continuity at {0} by 
setting U(0) = U(0) = (0,0)'. Moreover, if G has independent components, then W = 0. 

Proof of Proposition 6.1. Let us define the functions p n ^ by 

PnAn) = - ln(ff„(rW)) = - ln(F„((0, . . . , 0, n , 0, . . . , 0))), i = 1, . . . , d. 

The generalized inverse of p n ,i is given for < f < r n by 

P^i( f ) =int^r>0: 1 {X j , i >F i '-(i-T/r n )}>k„f^ =r n F i (X^ kn f]) ti ) 

since F*~ (F i (X^ knf ^ i )) = X( -i w . Without loss of generality, assume that p^~, (0) = 0. 
Note that p^(-) is a caglad function on [0,1]. Letting 7T;(t) =t<, we have 

ff n (-r) = ^(^(^(t)), . . . ,p^ d (7r d (r))). 

Let us introduce the functions p n , p™ v and from [0, l] d to M d defined by 

Pn(r) = (p„.i(7Ti(r)), . . . ,P n ,d(^d{T)))' , 

ed(r) = r. 

By Theorem 6.1, we have p n (-) =>• e<j(-) in D([0,l] d ) x •■ ■ x D([0, l] d ). It is easily deduced 
that p™(-) => e d {-) in D([0, l] d ) x ■ ■ • x D([0, l] d ). 
Let us define the processes 

V n (r) = y/k^(H n (r) - H{t)) 
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= V n (r) + ^/K (P(Nll\ = 0) - P(N^ = 0)), 
W n {r) = y/K{{-\nH n { T )) - (-lnff(r))) 

= W n (r) + \fK{r n {l - F(u rn (r))) + ln#(r)) 
and let U n (r) = (V n (r), W n (r))'. By (C3)(i), we have 

sup \^(P(N^ A = 0) - P(ivM =0))h0 



as n — ► oo , 

Te[o,i] d 

sup \\fk^{r n (\ — F(u rn (r))) + ln£f(r))| — * as fnoo 
re[o,i] d 

and it follows that 

in D([0, l] d ) x Z?([0, l] d ). By using the continuous mapping theorem (CMT) and similar 
arguments as in the beginning of the proof of Theorem 4.2 in [37], we deduce that 



in D([0,l] d ) x D([0,l} d ). 
Next, note that 



V n (r) = V n (p™(r)) + VK(H(pT(t)) H(t)), 
W n {r) = W n (pT(T)) + VM- lnH(p^(r)) (- mff(r))). 
Since W n =>W in D([0, l] d ), we have 

Vhr{Pn{-)-e d {-))^Z{-) 

in .D([0, l] d ) x • • ■ x D([0, l] d ). By using Vervaat's lemma [44], we get 

in D([0, l] d ) x • ■ ■ x D([0, l} d ). We deduce from the differentiability of H and H, and the 
finite increments formula, that 

H(tfT(.))-H(.) \ ( -VH(-yZ(-) 
lnH(p^(.))-(-lnH(.))J l^CJVW) 



in £>([0,l] d ) x L>([0,l] d ) 
Finally, we get 



[/„(•) ^ £/(•) + 



-wh(-)'z(.) 

H-\-)VH{-)'Z{-) 
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in D{[0,l] d ) x D([Q 1 l] d ). □ 

Let k = l/sup{(-lnff(x))- 1 \ff =1 n:T e [0,l] d \{0}} and introduce, for (k,t) G 
[0,x] x [0,l] d , 



W,f(x J T) = v^(Q!(x J T)-x) ! 



where 



ff n( x ' T ) = rE 1 W ( '" T) =0t and <>>' r ) = T~Tl N r 

3 = 1 3=1 



We define an additional intermediate -D([0, it]) x D([0, iir])-valued process by U^(>c, r) = 
(V£ (x,t), Wjf (x,r))', x: G [0,x]. Observe that t/f depends on the estimated series 
(Z^ T '')/>i and on the unknown threshold Vn It is worth mentioning that t/f and 

£/„ ar e closely related since 

U z n (x,T) = UM-lnH(T))- 1 r). 
Corollary 6.1. Suppose that (CI), (C2) and (C3) hold. Then for r G [0,l] d \{0} ; 

U*((-),T)=>tj(T{-]nH{T))-H-)) 

m £>([0,x]) x Z>([0,x]). 
Proof. We have 



isl- 
and it follows that 



V*i - iu("(-i»flW)" 1 T) 

^ {Ui=l d (^.-> X <r. 1 x<-.nH M ,-l fc „l)„)} JVr »J 



t7^(^r) = C/„(r(-lni7(x))- 1 ^). 



Fix r G [O,l] rf \{0} and consider the function xh [/^"(x, t) from [0,x] to M 2 as an 
element of D([0, X]) x D([0,x]). Since the map from D{[0,l} d ) to D([0,k}) taking /(•) to 
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/((-lnF(T)) _1 r(-)) is continuous for any r € [0, l] d \{0}, we deduce, by the CMT, that 

f/f((-),r)^J7(r(-lniJ(x)r 1 (.)) 
inO([0,x])xD([0,i<]). □ 

We now derive, from Proposition 6.1 and Corollary 6.1, the distributional asymptotics 
of the estimators. Let us define, for r £ [0, l] d \{0}, 

9 n (r) = ~ lDJ l" (T) , 6(t) = -4— (V-WIT-Ht) + W(r)8(r)). 
-\nH n (r) lntf(r) 

Note that O(-) has continuous sample paths on [0, l] d \{0}. 

By proper cube we mean a cube included in [0, l] d which does not contain {0}. 

Corollary 6.2. Suppose that (CI), ( C2 ) and ( C3 ) hold. Let Cd be a proper cube. Then, 
inD(C d ). 

Proof. We first recall that a map $ between topological vector spaces Bi, i = 1,2, is 
called Hadamard differ entiable tangentially to some subset S C B\ at / e £?i if there 
exists a continuous linear map V$(/) from £?i to B^ such that 

*(/ + ^)-*(/)_»yg (/) , g 

for all sequences t n | and g„ G £?i converging to g G 5. Let D{Cd, E) (resp., C(Cd, E)) 
be the space of functions from Cd to the set £cl which are "continuous from be- 
low, with limits from above" (resp., continuous). Let us consider the map <f> from 
D(C d , (0, 1)) x D(C d , (0, oo)) to D(C d , (0, oo)) defined by 

*(/i,/ a ) = ^. 

J2 

Note that this map is Hadamard differcntiable tangentially to C(Cd, R) x C(Cd,K) at any 
(fi,f-2) € C(C d , (0,1)) x C(C d> (0,oo)). Moreover, V$(/i,/ 2 ) is defined and continuous 
on C(C d ,R) x C(C d ,R) and is given by 

1 111 f 

V$(/i,/ 2 ) ■ (31,52) = --r-rgi + 

Jl/2 \J2) 

Since 
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we deduce by the 5-method (see Theorem 3.9.4 in [43]) and Proposition 6.1 that 

in D(C d ). □ 

We end this section with the proof of Theorem 4.1. 

Proof of Theorem 4.1. Let m > 1 and ti, — ,T m S There exists d C [0, l] d such 
that t\/L(ti), . . . ,T m / L(T m ) £ Cd- By Corollary 6.2, we deduce that 



^(^(Ti) - ^)) 4=1 ,..., m =>■ ( e (j^))) 



i— l,...,m 



Let m > 1 and r i , . . . , r m € ^fjc. By using similar arguments as for the proof of Corol- 
lary 6.1, we have 



(r^O.i-ow m => c; ' (•) 

.(- In H (Ti)) 



i—1 ,m 



in (D([0, k})) . Let us consider the thresholds 

v^(k) = Z^ - v (T<)( r n(-^H(Tj)) \ 

Recall that 



i = 1, . . . ,m. 



/ fc n (g^((.),r. t )-(-))^^ (■) 



(-lnF(r 4 )) 

in D([0, vc\). By using Vervaat's lemma [44], we deduce that 

^ ^ -(■))*-»(h " 



in £>([0, >r]). Note that 



3=1 



-9(Ti)(0 



'(rM-)i) 
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We then deduce by the CMT and Corollary 6.1 that 



V'(-lnH(T,)) 



in D([0, n]), by the finite increments formula and the CMT that 

v ^ (e -e(^)rn(-InH(x 4 ))/^ (01) _ e -9(T«)(0) e(Ti)e- e ^)(')^f (•) 5 ) 

V (-lnF(ri))/ 

and by the 5-method that 



Since 



-(•) V V W (-lnff(T0)7 ^ " V'i-lnHiT,)))) 

= J_( e o(^)()v( (.) 2 ) +6(Ti)w( (•) 

-(■)V V W (-lnF(ri))y V (-ln^r,)), 



lntffx =-x and H[x 5 = e" 9 ^, 

V (-lni/(r. t ))y V (-lnfr(TO)y 

we have 

— (e e ^v( h 5 ^ + 0(t;)w( x 5 ^ =e( x T 4 

Finally, fix k and deduce that 



7. Discussion 

In this paper, we have developed new estimators for the multivariate extremal index 
function. In order to construct scale invariant estimators, we have used a homogeneous 
transformation for the first estimator, but it leads to the question of the choice of the 
optimal transformation. We have also considered a second estimator which is scale in- 
variant without transformation. One may also exploit averaging methods and consider, 
for example, the estimator defined by 

Of (r) =lf* -^(to) dX) T e Uo s {XTQ : K > 0}) 
-lntf n (xT ) 
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where < a < <fi < oo. We have studied the weak convergence of our estimators as point- 
wise estimators and given their asymptotic distributions. To study their weak convergence 
as functional estimators, one must construct a specific functional space which is differ- 
ent from the Skorohod space of caglad functions which does not contain the set of scale 
invariant functions, then study their asymptotic properties in this space. This seems to 
be an important avenue for future research. 

Appendix A. Proofs of Propositions 2.1 and 4.1 



Proof of Proposition 2.1. Let A > 0. We have 



where (r) = inf {x E K : Fi (x) > t} . Since Ff~ (r) < FfZ. (r) for each r G (0, 1), we have 



nP(Z\ T) >n\- 1 ) 




= n 




0<nP U(*M>«»,iM)) )- nP [ \J{ x l,i>K-0--n~ 1 r i X)) 





< nP[ \J(F^(1 - n- l nX) < X lti <F£_(1- n-ViA)) 





\i=l 
d 



< 



J2nP(Xi, i = FIZ{l-n- 1 T i X)). 



Note that 



lim 



P(X lii = Fr-_(l-n- 1 n\)) 
nX/n 



= lim 



PjXj^F^l-n-^X)) 

P(x lii >u n! i(nX)) 

Fjjx) - Fj(x-) _ 
Fi(x) 



= lim 
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and then 

Urn nP(^ (x) > nA- 1 ) = Urn nP\\\(X u > u^faA)) ) . 

\i=l / 

By (1.2) and the homogeneity property of — lnlf, it follows that 

lim nP{z[ T) > n\~ l ) = lim nil - F{u n {\r))) = - \nH{\r) = A(- lnF(r)). 

ra — >oc n — >oo 

By taking A = k{— lni/(r)) _1 , we deduce that 

lim nP(Z ( , T) > v£\x)) = x. 

n — >-oo 

We now have 

pf max Z, W <«W(x) 

\ t= l,...,n 

= PI max max $ < n{— In H{t))x~ 

V /— 1 n i= 1 , d ' 



Note that 



= P\ max Tj max < n(— In if (t))x 1 

Vi=l,...,d Z— l,..,n ' 

= P( max -<n(-lnff(T)W-M 

v ' - 1 /•'.. (M n , t ) " / 

- P(M n ,i <F£_(1- n-^xni- \nH(r))- 1 ),i = l,...,d). 

P{M n ,i <F^_(\-n- X KTi{- lnP(r))- 1 ), i = 1, . . . ,d) 

- P(M n ,i < Ft {I - n^xnlr InH(r))- 1 ), i = 1, . . . , d) 
= P{M n<i =F£_(1- n-WiC- ln^r))- 1 ),? = 1, . . . , d) 

d 

< VnP(AT M = J F i ^(l-n- 1 xr i (-ln#(T))- 1 )) -> 0. 

It follows that 

lim P(M(*> <«W(x)) 

n— >oc 

= lim P(M„,i<Pr(l-n- 1 xn(-lnfl'(T))- 1 ),i = l,...,(i) 

n — »-oo 

= PT(>fT(-lnF(x))- 1 ) = F(^r(-lnPr(x))- 1 ) e ^ T (- ln ^( T » _1 ' 
= e -9(T)« 
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which means that 0{r) is the univariate extremal index of the stationary sequence 

(4 T) )n>l. □ 

Proof of Proposition 4.1. The proof follows the corresponding lines of the proof of 
Proposition 1 in [30]. Let Ti,T2 G [0,oo) d \{0} and s be a positive constant. Define 

f> („„ _ _ ) _ /A/-(Tl,T2) ^(1-1,1-2) »r(Tl,Tj) Ar(Ti,-T 2 )y 

By considering a A(m„(xi), ii„(T 2 ))-separating sequence, (r n )„>i, and Berstein's blocks 
method (see Lemma 2.2 in [15] or the proof of Lemma 6.7 in [37]), we get 

lim |E(e h/ - R "( sn '' ri >' r2 )) - (E(e il, ' ii " (r "' T1 ' 1 " 2) )) m "| = 0, 

n— »oc 

where m n = [sn/r n \ 5 L^J denotes the integer part of x and v € M 4 . Now, note that 

E(e iv ' R *^< T ^) =P(Nl2r n 2) =0) +E(e i ^" (r "^ 1 ^ 2) |^ X 2) >0)W^oX 2) >0 )- 
Since (r rl )„>i is a A(m„(ti), it n (T 2 ))-separating sequence and 

lim P{N { J^ 2) = 0) = c -e(-iV-2)in(i^ lW2 )) ; 

n— >oo ' ' 

we have 

lim — P(NZ£'r 3) > 0) = -6»(n V r a ) ln(i?(n V r 2 )). 

n^oo r n i i n 

We then deduce that 

E(e to ' fl »<"*' T1 ' Ta >) = exp(m, 1 P(7V^ 1 ^ 2) > 0)E(e to ' fl »< r »'' ri '' ra > - 1\N^;^ > 0)) + o(l). 
On one hand, we have 

lim E(e i "' jR " (sn '' ri '' r2) ) =exp(-«0(Ti Vr 2 )ln(F(ri V x 2 ))(Ee i^ '' c ^ 1 ' T2, - 1)) 

n — ► oc 

where = (Cg llTa) + Ci?'^ + Ci?" 2 \ C&' Ta) , C^Y- In particular, 

we derive the weak convergence of the sequence \Rn{n,T\,T2))n>i by choosing s = l. 
On the other hand, it is easily seen by using the definition of u n (r) that 



lim E(e h/i? '" (;in '' ri <' r2) ) = lim E(e iv ' R ° n( - sn > STl ' ST ^) = lim E(c k/jR ' l(n ' s ' ri ' s ' r2) ) 

n — ► oo n — >oo n—*oo 

= exp(- S 0(r 1 V r 2 ) ln(F(n V r^)^''^ 1 ^ - 1)). 
Therefore, Ee 1 " C (ST1,ST2) = 'Ee lv C (T1,T2> anc ] ft follows that n^ 1 ^ 2 ) is scale invariant. □ 
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Appendix B. Covariance function of 

Let us define the functions Ci.j(-, •) for i = 1,2 and j = 1,2 by 

<?!,!(-!-!, T 2 ) = H( Tl V T 2 ) - H( Tl )H(T 2 ), 

C 2>a (T ls T 2 ) = -6>(n v r 2 ) in(# ( Tl v r 2 ))E((C^^ 2) + C| Xl ' X2) )(d Tl,T2) + Cl Xl,X2) )), 

C 1 , 2 (r 1 ,x 2 )=E(^ 1 ^h {<1 ,. 2)=0i<1 ,. 2)=0} )+ J ff(x 1 )lniJ(x 2 ), 
C 2 , 1 (r 1 ,x 2 )=E(^ 1 ^h {<1 , T2)=0i<1 ,. 2)=0} )+ J ff(x 2 )lniJ(x 1 ). 

Note that Ci, 2 (t 2 ,ti) = C 2j i(ti,t 2 ). 

Let us now characterize the covariance function of the process 6. We have 

cov(0(Ti),e(T a )) 

Gi,i(ti,t 2 ) + - ^ T2 L r Gi,a(ri,r a ) 



H(Ti)Infl-(Ti)fl-(T 2 )Infl-(T 2 ) ff(Ti)Inff(Ti) 



H{T 2 )\nH(T 2 ) 



G2,i(i-i,t 2 ) + 6I(ti)6I(t 2 )G 2i2 (ti,t 2 ), 



where 



Gr.iCrr,^) = <?i,i(ti,t 2 ) - ( Tl )£ ^^C 2i1 (t« r 2 ) 

2—1 ' 



- Fi 2 

2 — 1 



+ if(r 1 )F(x 2 )^^^- Ca.a^ ,r a ), 

2—1 J — 1 



Gi. 2 (n,x 2 ) = d, 2 (Ti,T 2 ) - g(ri)E aing(Tl) C 2 , 2 (r< i;) ,T 2 ) 



d 



d d 



2—1 l — l ' 
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G 2 ,i(ti,t 2 ) = C 2 ,i(ri,r 2 ) + > — 5 C 2i i(ri ; ,x 2 ) 



ff ( T 2) aZ ^2, 2 (ri , r 2 ; ) 



9T ij2 

<2 d 



ff(T2) ^^^^ 9^ 2 ' 2( 1 ' 2 } ' 

t =i J= i ■ 



G 2 . 2 (ti,t 2 ) = C 2 . 2 (ti,t 2 ) + } — C 2 ,2V r l ,t 2 ) 



d\nH{r 2 ) (i) 

-02,2(Tl,T 2 J 



2 — 1 

y-V^ 91ni?(n) dln#(T 2 ) 



i=i j=i 



9r 



i, 2 



L-2, 2 (T 1 ,r 2 J. 
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